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Abstract 

Taking the standard zero curvature approach we derive an infinite set of inte- 
grable equations, which taken together form the negative Volterra hierarchy. The 
resulting equations turn out to be nonlocal, which is usual for the negative flows. 
However, in some cases the nonlocality can be eliminated. Studying the combined 
action of both positive (classical) and negative Volterra flows, i.e. considering the 
differential consequences of equations of the extended Volterra hierarchy, we deduce 
local equations which seem to be promising from the viewpoint of applications. The 
presented results give answers to some questions related to the classification of inte- 
grable differential-difference equations. We also obtain dark solitons of the negative 
Volterra hierarchy using an elementary approach. 

1 Introduction. 

The present paper is devoted to an integrable hierarchy which can be viewed as an ex- 
tension of the Volterra hierarchy (VH) by taking into account the so-called 'negative' 
flows. The idea of negative flows is not new. One of the first examples (and maybe the 
most striking one) is the relation between the AKNS and sine-Gordon models 0: the 
sine-Gordon equation can be derived as the negative flow of the AKNS hierarchy. Some 
recent results and approaches to this problem can be found in |21ISII11E1- The negative 
flows have been constructed for almost all known integrable systems. But there are a 
few exceptions, and the Volterra model seems to be among them: to our knowledge, the 
systematic derivation of its 'negative' equations has not been discussed in the literature. 

Some of the equations belonging to the extended Volterra hierarchy are already known. 
First, it contains the famous 2D Toda lattice ;6j. We would also like to mention the 
works of De Lillo and Konotop jZj , who discussed a nonlocal modification of the Volterra 
chain and recent papers by Boiti et al |Hl El, who derived and studied new integrable 
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discretization of KdV. As will be shown below, these equations belong to the extended 



Another indication of the fact that the neagtive Volterra flows should be studied more 
carefully is related to the problem of classification of the integrable discrete equations. In 
the paper by Yamilov ITU] it has been shown that all equations of the form 



where p, q and r are polynomials in n. This equation is known to satisfy integrability tests, 
to posses an infinite set of local conservation laws and generalized symmetries. However, 
as was stated, for example, in "it is the only example of nonlinear chain of the form 
([T]) which cannot be reduced to the Toda or Volterra equations by Miura transformations" . 
In this paper we will show that there is no need to look for substitutions converting (j2)) 
into Volterra chain, equations ^ (at least the simplest of them) are nothing but equations 
of the extended VH. 

The most straightforward approach to the negative flows is to consider them in the 
framework of the inverse scattering transform. In section |21 we start with the standard 
zero-curvature representation (ZCR) for the VH and derive equations of the negative 
Volterra subhierarchy. The resulting equations turn out to be nonlocal. It is usual for 
the case of negative flows. However, nonlocality does not mean that equations are not 
interesting from the viewpoint of possible applications. Moreover, sometimes nonlocality 
can be eliminated. This can be achieved by a redefinition of the variables or by considering 
some differential consequences of the equations in question. Namely this is the topic of 
section 121 where we present the local equations which in this or another way are related to 
the nonlocal Volterra equations. Finally, in section |3] we derive the soliton solutions for 
the negative VH. 

2 Zero curvature representation. 

The inverse scattering approach for integrable systems is based on the ZCR, when our 
nonlinear equations are presented as a compatibility condition for some linear system. 
For example, the Volterra chain. 



VH. 



can be divided into three classes. One of them can be presented as 



(1) 



(2) 



Un+1 ~ Un-l 



n = 0,±1,±2,... 



(3) 



is the compatibility condition. 




(4) 



for the system [T2] 
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Un^n 



(5) 



where 



c/„=t/.(A) = (_\ ;»), v.^v.(x)^{:i j;:^,) (6) 

and A is an arbitrary constant. Traditionally Volterra and Toda chains are considered 
more often in the framework of the 'big' Lax representation, when the system is finite 
{n = 1, N) and U and V are N x N matrices. In this paper we will use the 2 x 2 U-V 
pair (jni). This approach is equivalent to the N x N representation and sometimes even 
more convenient (that is, it can be more easily modified to the cases of different boundary 
conditions, including the soliton case when N = oo). 

The choice of the U-V pair for the system (jS)), even for a given ?7-matrix, is not 
unique. The matrix \^ in © is a second-order polynomial in A. By simple algebra one 
can find other matrices V, which are polynomials of higher order, such that ZCR (j3)) 
will be satisfied for all A provided u„ solve some nonlinear evolutionary equations. These 
equations are called 'higher Volterra equations'. Taken together they constitute the VH. 

Now we are going to derive the negative Volterra flows. The idea is simple: let us 
search for the V^-matrices which are polynomials not in A but in A~^. For example, the 
simplest of such matrices, which is linear in A~^, is given by 

Vn=(j_^] (7) 

\ Ap„_l Pn-l / 

It is easy to verify that (0) will be consistent if m„ and Pn satisfy the system 

= Pn-l Pn /g\ 
= Pn~lPn 

This is the first negative Volterra equation. 

Now our aim is to obtain an infinite set of the similar V^-matrices which solve ZCR 
(jH). Using the notation 



v;= j!" ) (9) 



r d 



we can rewrite (jD) as 



and 



= A (a„+i - a„) - bn+i - u^Cn, (10) 

= bn + UnCn+l, (U) 

= Xcn+1 + a„ - dn+i (12) 



tin = Un (a„+i - dn) - A6„ (13) 
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or, after eliminating 6„ and c?„, 

= A [ttn+l - an) + Un+lCn+2 " «nC„, (14) 
Un = Un [an+l - On-l + A (c„+i - C^)] (15) 

It is easy to show that (fT^ and (fT3j) possess solutions where are polynomials of the 
(2j — 2)th order in A^^ while Cn are polynomials of the (2j — l)th order for j = 1, 2, .... 
In what follows we indicate different polynomials with the upper index, al^\c!^\ and 
introduce an infinite set of times, t, , to distinguish the resulting nonlinear equations. 

By simple algebra one can establish the following relations between different polyno- 
mials: 

a(^+i) = X-'al^^ + X-'ali\ (16) 
= A-^c(f) + A-Si^-) (17) 

where a^^^ and 7^-'^ do not depend on A. Substituting (fTH|) and (fT7|) into (fT^ and (fT^ 

one can convert our system into 

a% - a(f ) + u^+ili^l - «n7i^'^ = 0, (18) 
7i^-) + a[f + = (19) 



and 



J-lnu„ = aJ2i-ai'|i (20) 



After introducing the tau-functions of the VH, 



Un = (21) 



equations (fTH|) and (fT^ can be solved. 



ai^) = ^In^ (22) 

Otj Tn 

1^ = T^ln-^ (23) 



dtj+i Tn-2 



and the resulting equations become 



Tn-iTn+i — In + - - lu r„ = (24) 



together with the first one 



Tn^iTn+i In - — = r„. (25) 

oil Tn-l 



These equations can be rewritten in the Hirota bilinear form as 
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and 



2Dj+i Tn+l ■ Tn-l + DjDi Tn ■ Tn = 



(26) 



Dl Tn+l ■ Tn-l — 



(27) 



where 



d - 

Dja-b=—a{...,tj + ^,...)b{...,tj-^,...) 



(28) 



5=0 



3 Some integrable models. 

In the previous section we have obtained an infinite set of integrable equations. These 
equations are nonlocal and it is difficult to say whether they are of any importance, say, 
from the viewpoint of applications. However, a hierarchy is more than a set of equations. 
The key point is that all equations of a hierarchy are compatible (some notes on this 
question can be found in the appendix). Hence we can consider them simultaneously and 
study the combined action of different ffows. In other words, we can take some finite set 
of equations of a hierarchy and analyse equations which follow from this system. So, this 
is the aim of this section. Starting with the equations of the VH (both negative, derived 
above, and positive, i.e. classical Volterra equations) we will deduce some their conse- 
quences, which are local and some of which seem to be more promising for applications. 



3.1 2D Toda lattice and the VH. 

First let us consider the simplest negative equation (|25|l together with the Volterra equa- 
tion (0) 

= - (29) 

d 

— lnM„ = Un+l-Un-l (30) 
Ot 



where p„ is related to r„ by 



Pn = 5 • (31) 



One can derive from the second equation that 

d 1 



which leads to 



r. Pn-i - Pn+l (32) 

UL Pn 



9^ 1 Tn-lTn+l - , 

^ lnr„ = — (33) 



dt dt Tn 
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where an unessential constant has been omitted. This equation, which can be rewritten 
in terms of p„ as 

lnp„ = -Pn-l + 2p„ - Pn+1 (34) 



dtdt 

is nothing but the 2D Toda lattice lE^. So, we have shown explicitly that tau-function 
()21|) of the VH is at the same time the tau-function of the 2D Toda lattice, or, in other 
words, that the 2D Toda equation can be embedded into the VH. The reverse statement 
seems not to be true: not every solution of (j33|) should solve (f29|) and (j30|) . In this sense 
splitting the (2+l)-dimensional equation into two (l+l)-dimensional equations from 
the VH is a kind of ansatz. But the class of solutions of the 2D Toda model which can 
be can be obtained by solving Volterra equations and (jHUj) contains a large number 
of important solutions such as, e.g., N-soliton and finite-gap quasiperiodic solutions. It 
is obvious that this is a relation between the 2D Toda model and the extended VH and 
could not be revealed if dealing with the classical Volterra equations only. 

3.2 Singular chain. 

Another consequence of is that the quantity 



d 

Xn = ^\'aTn (35) 

oil 



satisfies 



Xn = ^ (36) 

where the dot stands for the differentiating with respect to the first 'positive' Volterra 
time (i.e. time which corresponds to 0), t, Xn = dxn/dt. Thus we have come to an 
equation of the type mentioned in the introduction. So, equations probably can 
be embedded in the extended VH. Note that the quantity x„ cannot be expressed locally 
in terms of the 'positive' Volterra subhierarchy, without invoking the first negative fiow. 

Equation (jHUj) seems to be non-typical for physical aplications. However, it can be 
easily converted to a more usual form. Indeed, it follows from that the quantities 
and P„, given by 

Qn = o 1^ Pn = Xn+lXn (37) 

2 Xn-l 

satisfy the system 

'^Qn = p _V + P \ (38) 

2P„ = -coth(g„ + g„+i) -coth(g„ + g„_i) (39) 

which is a Hamiltonian system 
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with the standard Poisson bracket and the Hamiltonian given by 

^ = ^ E {111 I - I + In I sinh (g„ + g„_i) 1} . (41) 

3.3 Integrable dynamics of roots of a 3-polynomial. 

The following system is a 'toy' model which can be derived from ()36|) by imposing the 
periodic conditions Xn+2, = ^n- Noting that the product 

n = (xi - X2) {X2 - X3) (X3 - Xi) (42) 
is a constant and rescaling the time t ^ U ■ t one can rewrite ()36|) as 



(43) 



or 



Xi 


= [Xi - 


X2, 




- X3 


X2 


= {X2 - 




){X2 


- 3:3 




= {^3 - 




1(3:3 


- X2 


k = 


V'iXk) 




k = 


1,2, 



(44) 



where V is an arbitrary monic 3-polynomial, and its roots: 



3 

V{x) = l[{x-Xk). (45) 

k=l 

Of course, it is a very simple ordinary differential equation, which can be studied without 
using the inverse scattering transform machinery. In terms of the first symmetric function 
of the roots, a = Xi + X2 + X3, it can be rewritten as the stationary KdV, 

attt - Qaf = (46) 

(to make the following formulae more readable we indicate differentiatings with subscripts) 
or, after introducing a 'logarithmic potential' ijj, 

A/i^ = (47) 

as 

iJtttt4^-4:^m^t + 3tPl = (48) 
which is the expanded form of the simple bilinear equation 

l4 



Dfij-ilj = 0. (49) 
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3.4 3+1 dimensional example. 

The last example we want to discuss differs from those presented above. All previous 
equations considered in this sectioin were difference-differential systems. Now we derive 
from the VH some partial derivative equations in (3-|-l)-dimensional (physical) space-time. 
Out starting point is equation which can be presented as 



Pn — In = — - - inTn- (50) 

dtj+i Tn-1 dtjdti 

Acting on both sides of this equation by the d/dik operator one can note that right-hand 
side is symmetrical in j and fc, which leads to 

( ^ ^ ^ ^ \\ ^"+1 n 

Pn 7T= 7^ Pn = U. (51) 

\dt.j dtk+l dtk Otj+iJ Tn-l 

On the other hand, differentiating ()50|) with respect to the first 'positive' Volterra time t 
one gets 

which can be rewritten using as 

at Otj+i Tn-l Otj 

It is easy to note that equations (jKT|l and form a closed system 

d dA d dA ^ 

^ P T^P 7^ = 

otj otk+i 01 k Ctj+1 

d OA dp 

P 7^ = 7^ (55) 



dt dtj^i dtj 



for two functions 



p = Pn and A = In -^iil n = constant (56) 

with n being fixed. 

From this multidimensional system one can derive a rather interesting consequence in 
3+1 dimensions. To this end we denote 

(x, y, z) = (ti, t2, h) (57) 

and introduce vector u by 

^ fdA OA dAY 
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In these terms equation ()55|) can be rewritten as 



Vp 



(59) 



where V is the three-dimensional gradient, V = {d/dx,d/dy,d/dz) , while equation ()54|1 
gives 

curl pu = 0. (60) 
Eliminating p from (j59|l and (j60|l one obtains that vector u satisfies the following equation: 



curl u 



U X 



du 
'dt 



(61) 



where [... x ...] is the usual three-dimensional wedge (vector) product. 

Equation (jHH), which can be termed as 'dual to' the Euler equation, du/ dt = [u x curl u] 
has already been discussed in the literature. For example, it has been proposed as a 
relativistic generalization of the curl Vg = condition for the velocity of a superfluid con- 
densate (see, e.g., [121 )• However, its relation with the integrable Volterra system (which 
also implies integrability of this model) seems to be new, and we are going to present 
more elaborated studies of in future papers. 



4 Soliton solutions. 

Here we would like to recall that equations are not only to be derived or classified, but 
should also be solved. In this section we present some class of solutions of the negative 
VH, namely the dark soliton ones. 

There are many ways to derive pure-soliton solutions for an integrable model: the 
inverse problem, the dressing method, Backlund-Darboux tansforms, Hirota's ansatz. The 
key idea of the last approach is that soliton solutions of all the integrable models possess 
the same structure, and the only thing one has to do while solving a particular equation 
is to determine some constants. In this paper we also exploit this fact. However, we do 
not use Hirota's technique but solve our equations using the elementary matrix calculus. 
Our starting point is that usually A^-soliton solutions can be constructed of determinants 
of some combinations of the N x N matrices having the following structure: 

(An),, = (62) 
i.e. matrices which satisfy the equation 

LAr,-AnR=\i){an\ (63) 
Here we use the 'bra- kef notation, 

(a„| = (a„,i, ...,a„,Ar) \ £) = {ii, i^) (64) 
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and L and R are some diagonal matrices 

L = diag (Li, L^) R = diag Rn) 



(65) 



4.1 Algebra of matrices ( 1621 ). 

Before proceeding further we present some formulae related to matrices ()62|) which we 
need to construct soliton solutions of the negative Volterra equations under the so-called 
'finite-density' boundary conditions. We do not consider the most general case of 
but restrict ourselves to 

L = R-^ (66) 
and specify from the beginning the n-dependence by 

A„+i = Ar^R (67) 

(as will be shown below, namely these matrices we need for our purposes). From (|63|1 . 
which we rewrite now as 

R-'An-An+i = \i){an\ (68) 

or 

A^Z-ZA^+, = Z\i){an\Z (69) 

where 

Z=(/ + i?-i)"^ (70) 
and I is N X N unit matrix, one can derive that matrices Fn inverse to I + An, 

Fn = {I + A„)-' (71) 

satisfy the identities 

Fn+iR~^ — R~^Fn-i = F„+i| £)(a„_i|F„_i (72) 

and 

- FnZ = ZFr,+i\£){an\FnZ. (73) 
These relations, together with the following formulae 

l-{an\FnZ\i) = ^ (74) 
l + (a„_i|ZF„|£) = ^ (75) 

where 
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ujn = det (/ + An) (76) 

lead to 

u^F^Zli) = LUn+^ZF^+,\i) (77) 

{an\ZFn+l. (78) 

Formulae ()74j) and (ffH|l follow from (jHTj) and the fact that for any bra-vector (i.e. N-tow) 
{u\ and any ket-vector (i.e. A^-column) |t>), det (/ + \v){u\ ) = 1 + {u\v), while fjTTjl and 
(ffHjl can be obtained from (ffHjl by multiplying it by \i) and (a„|. Formulae (ff^ - (ffHjl are 
all we need in the following consideration. 

4.2 Auxiliary sytem. 

Now we are ready to solve auxiliary equations 

8,11,"^ = —^ 1 (79) 

and 

2 

dj+i In ^ = ^ dij In tOn. (80) 

Solutions of these equations then can be easily modified to become solutions of the negative 
VH under the non-vanishing boundary conditions — > ±oo 

The crucial point of the ansatz (j62j) is that dependence on times ii, i2, ... is incor- 
porated in the A^-rows (a„| = (a„ (ti, •••) |- Moreover, we assume (which is usual for 
pure soliton solutions) that In a„fc is a linear function of times. Hence, differentiating the 
A- matrices leads to multiplication from the right by some constant matrices, Cj, 

^An = A^C, (81) 
which should be determined. As follows from (jHT|) . 

^ hiujn = tr F^A^Cj = tr (/ - F„) C, (82) 



which gives 



|_ In = tr (F„_i - F„+i) C,. (83) 



Using (f?^ and the fact that R and Cj commute, the last equation can be rewritten as 

d UJn+l 
dtj U!n-1 



^ Ir,'^ = tr - C,i? (84) 



-K-i|i^n-iC,i?F„+i|£) (85) 
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On the other hand, by multiplying ()74|) and ()75|) one can derive the identity 

1- (a„_i|F„_iX„F„+i|£) (86) 



where 

X„ = Z (/ + A„+i) - (/ + A.-i) RZ + Z\i){a^\Z (87) 
Applying (f72j) one gets 

Xn = (I -R)Z = constant (88) 

i.e. 

,2 



- 1 = {an-i\Fn-iiR - I)ZFn+i\(.) (89) 
Comparing and the last formula one can conclude that Un solve (jZHI) if 

= - /) Z (90) 

or 

Ci = {I ~R){I + R)-\ (91) 
Now our aim is to solve (jHUjl . Differentiating (|H^ with respect to ti we get 

-^-^ In Un = tr (92) 

OtiOtj 

From definition ()91|) of Ci and ()69|) one can derive 

Ar,Ci = ZA^~Ar,Z + Z\i){an-i\Z (93) 

which leads to 

A„F„CiF„ = F„Z - ZFn + FnZ\i){a^_,\ZFn. (94) 

Noting that Z and C-,- commute and that the trace of a commutator is zero, the right-hand 
side of fl92j) can be rewritten as 

tr A„F„CiF„C,- = {a^.,\ZF^CjF^Z\£) (95) 
which gives, together with (fTTjl and (fTHjl . 

ln,;„ = ^"-^t'"+^ (a„._i|F„_iQZ^F„+i|£). (96) 

Comparing this expression and equation (|H^ for the derivative with respect to the (j-l-l)th 
time, ij+i, one can conclude that Un is a solution of (jHUj) . if matrices Cj+i and are 
related by 
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Cj+1 = CjZ^R"\ (97) 
This recurence together with the 'initial condition' (jnT|) can be easily solved: 

In such a way we have shown that if the time-dependence of the A-matrices is given 
by (|RT|l and (jHE)), then the quantities Un given by (|7H|l solve (|Hn| . 



4.3 Dark-soliton solutions of the VH. 

It is easy to see that equations (f?^ and (jHUj) differ only in a constant term on the right- 
hand side of (f7n|l , which can be taken into account by exp (nti/2) multiplier: functions 



r„ = exp 



nti 



LOr. 



(99) 



solve ()25|) . The limiting values of the coresponding u-functions are unity and to meet 
general 'finite-density' boundary conditions, 



lim Ur, 

n— >itoo 



(100) 



one has to add the u^^'^ factor and to rescale the times, tj tj/ul^ The final formulae 
for the iV-soliton solutions of the VH can be written as 



Tn = U 



n74 



nil 



exp 



det 



-L -r — I 



Rk 



(101) 



j,k=l,...,N 



The time dependence of on negative times is given by 



ak{t) = ak (ti, ^2, •••) = af^ exp ( ^ u^f' tj 



(102) 



with 



1 1-Rk 

^Uoo 1 + Rk 



(2 + Rk + Rk ^ 



(103) 



where Rk, a^k^ , k = 1,...,N are arbitrary constants. Note that constants ij appearing 
in (j62|) have been incorporated in af^ by the transform An M'^A^M with M = 
diag (£1, ...,iN) which does not change determinants. 
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5 Conclusion. 



In the present paper we studied the negative Volterra flows. The main result of this work 
is that even in the theory of, so to say, classical systems there are some questions which 
have not been discussed in the literature yet. We applied the widely used zero curvature 
approach to the well-known scattering problem and obtained some results which seem 
to be new. Deriving equations and or and (j77j) is not a very difficult 
problem, though it fills some gaps in the theory of one of the oldest integrable systems. 
However, a hierarchy is more than a set of equations, and one can hardly enumerate 
all equations which are contained in this or that hierarchy, i.e. to study all possible 
differential consequences of equations of a given hierarchy. In section |2] we demonstrated 
that from the extended VH, which is a set of differential-difference equations, some of 
which are nonlocal, one can extract rather different systems. For example, in our opinion 
it is very difficult to say a priori that equation (a vector partial differential equation 
in 3+1 dimensions) has any relation to the Volterra chain. In future papers we are going 
to present some other integrable systems which can be 'embedded' into (or extracted 
from) the extended VH, and we hope that the results presented above demonstrate that 
the theory of integrable systems is far from finished. 
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Appendix. 

The question of compatibility of the equations of a hierarchy (or, in other words, of com- 
mutativity of the corresponding flows) is of vital importance and usually not trivial. The 
most straightforward approach to this problem is to derive the Hamiltonians generating 
these flows and to show that they are in involution. However, the Hamiltonian repre- 
sentation of the negative Volterra hierarchy is still unknown and is a subject of future 
studies. Thus here we do not discuss this rather serious problem in the most general way 
and restrict ourselves to proving the commutativity of the first (classical) Volterra flow 
Q and all the negative flows derived in this paper, a fact which is crucial for what has 
been presented in section |3] 

This claim can be verified directly in the simlest cases. Let us first consider the 
Volterra equation (jSI), rewriting it now as 



which in its turn leads to 



J-ln^ = «„, (104) 

dti Tn-l 



In = Un+l + Un, (105) 

oh Tn-l 
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and show that it is compatible with 

^ln^ = -. (106) 

dti Tn-l Pn 

To do this we need auxihary formulae for the derivatives of Un and p„. The first one 
(which follows from ()106|) ) is 

= fln^-l„^)=l-^ (107) 

Oh dti \ Tn-l r„_2/ Pn Pn-1 

while the second one is a consequence of ()104|1 : 

7^ lnp„ = 7^ (in^^ - In^^ I = - (108) 

dti Oti \ Tn Tn-l) 

Now it is easy to derive 

d ( d ^ Tn+l\ / 1 1 \ / 1 1 \ , , 

^ — In = Un+l \+Un\ = -Pn+l + Pn^l (109) 

dti V^^l '^n-lj \Pn+l Pn) \Pn Pn-1 J 

(we have used p07|) and the fact that m„ = PnPn~i)- 
On the other hand, it follows from ()106p and ()108|l 

( i In ) = — (m„ - Un+l) = Pn-1 " Pn+l- (HO) 

dti \dti Tn-l J Pn 

Since the right-hand sides of p09|l and piOj) are equal, one can conclude that 

dd^^r^^dd^^Tn^ (111) 

dti dti Tn-l dti dti Tn-l 

which proves the compatibility of (|104p and (jl06|) . 

Further we proceed by induction. Suppose that Volterra flow ()104|) commutes with 
all the negative flows d/dii, d/di2, ... , d/dij. Now our aim is to show that it commutes 
with the (j + l)th negative flow as well. Indeed, differentiating equation (j21I), which can 
be rewritten as 

^ ln^ = -l-i^lnr„ (112) 



dtj+i Tn-l Pn dtidtj 

with respect to ti, one can obtain 



dti \dtj+i Tn-l J dti \pn dtidtj / 

-(——]-^\nT-——(-^\nT 

\dti Pn) dtidtj Pndtj \dtidii 
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Id 

= {Pn+l - Pn-l) r.T + — (113) 

OtiOtj Pn Otj 

Here we have used the commutativity of d/dti and d/dtj as well as the identity 

32 



which can be derived by differentiating (|1U6|) with respect to ti (see section . Noting 
that 

d d 1 d"^ T 

^Pn = -pl^— = -pI t^tw' ^ 

otj otj Pn OtiOtj Tn-1 

one can rewrite pi3p as 

-K- T^f — In^^ = {pn+i -Pn-i) In + In (116) 

dti \dtj+i Tn-iJ OtiOtj dtidtj r„+i 

On the other hand, 



d f d , Tn+l\ d 



77- In = T^r {Un+l+Un). (117) 



^j+l Tn-l) dtj+i 

One can obtain from (|112j) the following formula for the derivative of 



Un = PnPn-lT^ ^ lu ■ 



rvi fl' IT- HIT' li — _L r\, 1 I 

= a;r I'^^n-i ~ ar a;r ^^^^ (^^^ 
atiOTj at 1 Otj 



Substituting this expression in ()117p leads to 



In = [Pn+l - Pn-l) rvr + P« -vr oX (l^^) 



OTj+i yOTi Tn-i J OtiOtj OtiOtj Tn+l 

Comparing the right-hand sides of pi6|) and ()119|) one can easily see that they coincide, 
which means that 

Oti dtj+i Tn-i dtj+i dti Tn-l 

This leads by induction to the fact that all the negative Volterra flows d/dij commute 
with the classical one, 8/ dti. 
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